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Abstract: Results from the new Monte Carlo event generator Herwig++ are pre- 
sented. This first version simulates Hadron Emission Reactions With Interfering 
Gluons in electron-positron annihilation. The parton shower evolution is carried out 
using new evolution variables suited to describing radiation from heavy quarks as 
well as light partons. The partonic final state is fragmented into hadrons by means 
of an improved cluster hadronization model. The results are compared with a wide 
variety of data from LEP and SLC 
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1. Introduction 

The new generation of liigli energy colliders such as the Large Hadron Collider (LHC) 
or a future linear collider (NLC) require new tools for the simulation of signals 
and backgrounds. The widely used event generators HERWIG [IH and PYTHIA 
underwent tremendous development during the LEP era and have reached the limit 
of reasonable maintenance in the future. Therefore these programs (Pythia?) as 
well as new projects, like SHERPA are being completely (re-) developed in the 
object-oriented programming language C++. 

In this paper we present results from the new Monte Carlo event generator 
Herwig++ as the first step in the redevelopment of HERWIG. The generator will be 
used here to simulate e^e~ annihilation events. In order to have full control of the 
basic physics steps that are simulated, we need to put the new generator on a firm 
basis with respect to LEP and SLC results before we go on to upgrade it to initial- 
state showers and the other requirements for the simulation of lepton-hadron and 
hadron-hadron collisions. Therefore we have tested the predictions of the generator 
against a wide range of observables that have been measured at LEP and SLC, and 
have explored the sensitivity to the most important parameters and cutoffs. We did 
not perform a high-precision tuning: our aim here is rather to describe the program 
and to show that it is able to give results as acceptable as those generated by its 
predecessor HERWIG for a reasonable choice of parameters. 

2. Main features of the code 

The details of Herwig++ will be described in conjunction with the release of the 
code 0. The main stages of the simulation of e^e~ annihilation are the same as in 
HERWIG However, in comparison to its predecessor, Herwig++ features a new 
parton shower and an improved cluster hadronization model. At present, hadronic 
decays are implemented in the same fashion as they were in HERWIG. 

The program is based on the Toolkit for High Energy Physics Event Generation 
(ThePEG) and the Class Library for High Energy Physics (CLHEP) [§. They are 
utilized in order to take advantage of the extended general functionality they can 
provide. The usage of ThePEG unifies the event generation framework with that of 
Pythia?. This will provide benefits for the user, as the user interface, event storage 
etc. will appear to be the same. The implementations of the physics models, however, 
are completely different and independent from each other. 
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Our simulation starts with an initial hard process e+e^ (7*, Z^) —^qq + 77. 
The final state photons simulate QED radiation from the initial state, so that a 
radiative return can be properly simulated. For the present paper we shall only be 
interested in the details of the QCD parton shower in the final state. The final-state 
parton shower starts with a quark and antiquark that carry momenta Pq and Pq, 
respectively, and have an invariant mass squared of = {pq +^5)^. The only detail 
we are concerned with in relation to initial-state radiation is that the centre-of-mass 
frame of the gg-pair is slightly boosted with respect to the collider laboratory frame 
and that Q may be different from the e~^e~ centre-of-mass energy. We have made 
sure that the applied cuts on the energy of the annihilating e'^e~ subsystem are the 
same as those used in the experimental analyses. 

2.1 Parton shower 

The partonic evolution from the large scale of the hard collision process down to 
hadronic scales via the coherent emission of partons, mainly gluons, is simulated 
on the basis of the Sudakov form factor. Starting from the hard process scale Qo, 
subsequent emissions at scales Qi and momentum fractions Zi are randomly generated 
as a Markov chain on the basis of the soft and collinear approximation to partonic 
matrix elements. Details are described in chapter 5 of ||^. In Herwig++ we have 
chosen a new framework of variables, generically called (g, z). Here, g is a scale that 
appears naturally in the collinear approximation of massive partonic matrix elements 
and generalizes the evolution variable of HERWIG to the evolution of massive quarks, 
z is a relative momentum fraction; the evolution is carried out in terms of the Sudakov 
decomposition of momenta in the frame where the respective colour partners are 
back-to-back. As in HERWIG, the use of the new variables allows for an inherent 
angular ordering of the parton cascade, which simulates coherence effects in soft gluon 
emission. The details of the underlying formalism are described elsewhere . 



The most important parameter of the parton shower that we will be concerned 
with in this paper is the cutoff parameter Qg, which regularizes the soft gluon sin- 
gularity in the splitting functions and determines the termination of the parton 
shower. Less important but relevant in extreme cases is the treatment of the strong 
coupling constant at low scales. We have parametrized as{Q) below a small scale 
Qmin > ^QCD in different ways. We keep Qmin generally to be of the order of 1 GeV, 
where we expect non-perturbative effects to become relevant. Below that scale as{Q) 
can optionally be 

• set to zero, as{Q < Qmin) = 0, 

• frozen, asiQ < Qmin) = "5 (Qmin) ) 

• hnearly interpolated in Q, between and a5(Qmin), 

• quadratically interpolated in Q, between and ^^(Qmin) ■ 
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We put the final partons of tfie sfiower evolution on their constituent mass shells, 
since the non-perturbative cluster hadronization will take over at this scale, so we 
usually have kinematical constraints that keep Q above Qmim in which case the 
treatment below Qmin is irrelevant. Typically, OisiQrain) ~ 1 here. 



2.2 Hadronization and decay 

The partonic final state is turned into a hadronic final state within the framework of 



the cluster hadronization model of HERWIG 11 . In order to address some short 



comings 



a new cluster hadronization model ^ has been created for Herwig++, 
which is discussed in sec. ^. The emerging hadrons are possibly unstable and eventu- 
ally decay. The decay matrix elements and modes correspond to those in HERWIG. 



3. The parton shower in detail 
3.1 Hard matrix element correction 

Before we begin the parton shower evolution, but after obtaining the final state qq- 
pair from the hard process, we decide whether or not a so-called hard matrix element 
correction will be applied. In order to do so, we decompose the gg^f-phase space into 
regions that will be covered by the parton shower emissions and a 'dead' region that, 
based on our choice of evolution variables and initial conditions, is never populated by 
first parton shower emissions (see ||10|)- To take into account gluon emissions into the 
dead region we generate a pair of three-body phase-space variables x, x according to 
the first order QCD matrix element. However, we only accept emissions into the dead 
region of phase space at a rate that is given by the QCD matrix element, that is, only 
3 % of all emissions are corrected by the hard matrix element at all. Once we accept 
an additional hard gluon emission, we replace the gg-final state with the qqg final 
state. We keep the orientation of either the quark or antiquark with weights and 
respectively, resulting in properly oriented three-jet events apart from finite mass 



effects [0]. In this way, we take into account the most important subleading higher- 
order corrections that are not enhanced by additional soft or coUinear logarithms. 



3.2 Initial conditions 

Having completed the hard matrix element correction, the next task is to determine 
the initial conditions for the parton shower evolution. For every particle a we de- 
termine the colour partner or, more generally, the gauge 'charge' partner a. In the 
case of a qq final state there is no ambiguity, whereas the gluon in qqg is assigned 
the quark or the antiquark with equal probability. 

new cluster hadronization model that addresses some of these shortcomings is also a feature 
of SHERPA [|13). 
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For different interactions there can be different 'charge' partners. In our case we 
have also implemented collinearly enhanced photon emission from charged particles. 
In the case of the qqg final state the gluon doesn't radiate photons and the only two 
charge partners are the quark and the antiquark. The remaining parts of the shower 
evolution are carried out in exactly the same way. Different sorts of interaction just 
add another splitting possibility for a given particle, which will compete with the 
others for the next possible splitting that occurs. 

Once the colour partners are determined, we fix the shower kinematics and the 



initial evolution scale. As explained in detail in |T^, the shower evolution of a particle 
a is carried out in a Sudakov basis, 

q = apa + + qj_, (3.1) 

where Pa is the momentum of particle a with (current) mass-squared p'^ = rri^, ria is 
a lightlike vector in the 'backwards' direction, along the momentum pa of the partner 
a, and g_L is the transverse momentum, q_L ■ Pa = q± ■ ng, = 0. In the centre-of-mass 
frame of a and a we have Pa = ^Q{l,v) and we set Ug, = \Q{v, —v), qj_ = {0,q±) 
with q± ■ V = 0. Given this basis, we calculate the initial evolution scale for each 
particle as 

Ci = (Pa + Pn) ■ (Pa + rig) = ^Q\l + v) . (3.2) 



We note that this is the most symmetric choice of initial conditions (see [0). In the 
gg-case, this choice starts the evolution of quark and antiquark at the same scale. 
We could as well choose another pair of evolution scales. If we do so, however, we 
make sure that the phase space region of soft gluon emission is covered uniquely 



and smoothly with the radiation from the two partners |T^. For later kinematic 
reconstruction we have to store the momenta pi of the outgoing partons at this 
stage. 

3.3 Parton splittings and kinematics 

Starting from the evolution scale (ji = gini we now carry out the parton shower 
evolution for each final state particle separately. For every possible splitting a ^ be 
of particle a we determine the scale of the next branching gj+i based on the Sudakov 
form factor 

Abaiq^Qi 



SbaiQi, Qi+l) = . - — T , (3.3) 

Aba{qc,qi+i) 



where 



^rfg^ f as{z,q) 



Abaiqc, g) = exp <; - / J dz A,(z, g)e(px > 0) !> . (3.4) 



qc is the lower cutoff of the parton shower which, by default, is taken to be the 
nonperturbative gluon mass = 750 Me V. as{z,q) is the running coupling in the 
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case of QCD evolution and generally depends on the evolution scale and momentum 
fraction. We choose ^(1 — z)q as the argument of the running coupling which re- 
duces to the transverse momentum in the massless case. Pba{z,q) are the quasi 
collinear splitting functions that depend on the evolution scale in the case of massive 
partons [|T5|. For QCD branchings they are 



Pqq{z,q) = Cf 



1 + 



2ml 



z{\ — z)q 



n2 



Pgg{z,q) = Tj 



R 



2z(l -z) + 



2ml 



zl 



z)q'^ 



Pgg{z,q) = Ca 



1 - z 



+ + z{l-z) 



Similarly, for the branching q 



q'j, ignoring the parton mass, we have, 
,l + z^ 



Pl{z,q) 



qq\ 



(3.5) 
(3.6) 
(3.7) 

(3.8) 



with Ca being the electric charge of the parton in units of the elementary charge. 
Of course we have to take the fine structure constant aem in eq. ( |3.4| ) in this case. 
Q{qi_ > 0) limits the phase space to the region where it is possible to reconstruct the 
transverse momentum from the evolution variables (g, z) , which is a complicated 
and implicit function in our case. However, using the veto algorithm described below 
we do not need to know the phase space boundary explicitly. 

The evolution variables q and z determine the kinematics of the parton shower. 
The momentum fraction z is simply the ratio of the Sudakov variables a in eq. ( |3.1j ) 
for the parent and daughter parton. 



zaj 



(3.9) 



Based on its meaning in the quasi-collinear limit, q determines the relative transverse 
momentum as 



\P±\ 
\P±\ 



\/z^(l — z^q"^ — jj? 



(quark branching) 
(gluon branching) 



(3.10) 
(3.11) 



where 

fi = max{ma,Qg) , (3.12) 

when a quark of mass a is involved in the branching and simply fi = Qg for the 
splitting g gg. Here we have introduced the cutoff Qg in order to regularize the 
soft gluon singularities in the splitting function. The relative transverse momentum 
p± is related to the Sudakov variables ( p.l|) of the parton branching as 



P± = q±i+i - zq±i 



(3.13) 



- 6 - 



From eqs. (|3.10 ) and we immediately get the phase space constraint for 



p± in eq. (|3.4| ). We require z to correspond to a real value of p±. For gluon splittings 
we explicitly obtain 



z- < z < z+, z± = ^yL±^l - ^ j and q > Afi . (3.14) 

For quark splittings the phase space boundary is the solution of a cubic equation but 
the allowed z range always lies within 

^<^<1_^. (3.15) 
q q 

Therefore it is simplest to generate z within this range and reject those values that 
lie outside phase space. Finally, it takes an azimuthal angle 0, which is currently 



chosen randomly and may later be related to spin correlations [IG], to complete the 



four-momenta of the parton shower in a final kinematic reconstruction. 
3.4 Parametrization of Qg 

We introduced the cutoff parameter Qg quite naturally as a gluon virtuality into 
the shower kinematics. Considering the phase space that is available to the parton 
shower, we would expect a natural threshold in q of the order of rriq + Qg. In 
contrast, we find from eq. ( p. 101) that the actual threshold behaves approximately as 



Qthr = 1.15(mg + 2Qg). Hence, we find that, particularly for heavy quarks the phase 
space limit is well above our expectation. 

There is no reason why Qg should be kept as the same parameter for all quark 
flavours. Therefore, we have chosen to parametrize the threshold for different flavours 
in terms of a unique parameter 6 as 

Q, = . (3.16) 

which leaves us with a threshold Qthr = 0.85 + 6 for all flavours^. The resulting 
phase space in q and z is then as shown in fig. 0. We show the dependence of 
our results on the parameter 5 in most of the plots of sec. ^. In the case of gluon 
splitting rriq is the mass of the splitting products, i.e. the quark mass in case of a 
g ^ qq splitting or rriq = in g gg splitting. 

3.5 Single branching process 

For timelike (i.e. final state) branchings, given an initial scale qt, the Sudakov form 
factor eq. (|3.3|) gives the probability for no branching above the scale qi+i- Hence, 
1 — Sab{(li,Qi+i) is the probability for the next branching to happen above g^+i and 



^In principle, the coefficient of m„ could be a model parameter. 
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Figure 1: Available phase space of light (left) and 6-quarks (right) for q ^ qg splitting 
for various values of Qg and depending on the parametrization in terms of 6, eq. ( 3.16| ). 
The dashed lines on the right correspond to the same Qg values as for the light quarks. 



its derivative with respect to g^+i is the probability density for the next branching 
to happen at the scale gj+i. 

We sample the next branching scale with the veto algorithm. We overestimate 
the integrands as follows. We take the absolute maximum of the relevant coupling 
«max as this is generally a very slowly varying function. The splitting functions are 
overestimated by 



2Cf 



9gg{z) = Ca 



I- Z 

Tr, 



1 1 
+ - 



z z 



2el 
1-z 



(3.17) 
(3.18) 
(3.19) 

(3.20) 



in such a way that their integrals Gba{z) are invertible functions. Furthermore, they 
do not depend on q anymore. The phase space in z is overestimated by taking the 
maximum value of the evolution scale, gj. From this, we calculate the limits in z 
from eq. ( |3.14| ) or ( |3.15| ), respectively. As we can now easily integrate and invert the 
exponent in eq. ( p.4| ), we can sample values and Zg- Then we subsequently apply 
vetoes with weights 



wi = Q{p± > 0) , W2 



Pbajz, q) 
9ha{z) 



a{z,q) 



(3.21) 
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When all vetoes are passed, we have a scale g^+i = Qs and a momentum fraction Zs 
value. If not, we try to obtain a new branching, now starting at scale g^, repeating 
until we either accepted a scale as the next branching scale or we obtain a scale 
Qs < Qc at which we cannot resolve a parton any further. 

In this way we calculate branching scales g^+i for every possible splitting process 
for a given particle. The splitting with the largest scale of those above Qc is then taken 
to be the next branching. In this way we can easily include any type of branching. 

3.6 Angular ordering 

Once a parton is split its resolution scale is still above the smallest resolution 
scale Qc- In order to have angular ordering we now calculate the subsequent branch- 
ings of its daughters as g^+i and fcj+i with the conditions 

Qi+i < ZQi ki+i < (1 - z)qi . (3.22) 

This branching process is repeated until no more daughter particles are resolved at 
scales above the resolution scale qc- Note that, for our choice of evolution variables, 
the parton shower is terminated because there is no more phase space available at 
low scales. The lower limit of evolution is normally given by the soft gluon cutoff 
Qg (or 6) or the masses m of the branching particles (cf. fig. |l|). However, when we 
choose very small cutoffs Qg, which are in principle allowed, we apply the additional 
constraint q > qc on the shower termination. 

3.7 Soft matrix element correction 



As explained in sec. we explicitly populate the 'dead region' of the qqg phase 
space according to the correct QCD matrix element. We also improve parton shower 
emissions within the shower regions of the phase space, as the parton shower might 
generate relatively hard gluon emissions which are not within the domain of validity 
of the quasi-collinear approximation anymore. 

In order to do so, we keep track of the relative transverse momentum p± (cf. 
eqs. ( p.lOp , (|3.11| )) that was generated during the parton shower evolution of one jet. 



Whenever we find that this transverse momentum is the largest that was generated 



during the evolution so far we apply a so-called soft matrix element correction W7 . 

We consider all previous gluon emissions as being infinitely soft in comparison to 
the one we are testing. This allows us to compute the three-body (i.e. qqg) variables 
(x, x) from the parton shower variables (g, z) and the respective Jacobian. Then we 
compare a random number with the ratio of the true matrix element to the parton 
shower approximation and reject the branching if the ratio is smaller. Clearly this 
requires the parton shower emission probability to be larger that the matrix element 
everywhere in phase space, which is true for this process and our choice of evolution 
variables. 
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3.8 Reconstruction of kinematics 

As we generate the parton splittings i {i+l)kwe can calculate the ttj+i-component 



and the transverse momentum of the daughter parton {i + 1) using eqs. ( |3.9|) and 



( |3.13|) . The respective variables of the second particle k are simply obtained using 



1 — z and the difference of the transverse momenta of i and (z + 1). However, we 
can only reconstruct the f3 variables when we know the virtuality of each particle. 
This is done recursively once the parton shower evolution has terminated. The final- 
state particles are put on their constituent mass shells and then we obtain the beta 
components from 

P = l±A^ , (3.23) 

Zap ■ n 

where p and n are the Sudakov basis vectors of the shower. These were determined 
in the initial phase and remain fixed for each jet. 

After the completion of the shower evolution of every parton j involved in the 
hard process, the jet parent partons are not on their mass shells = anymore. 
Instead, they have acquired virtualities g|. If the original momenta were given as 

Pj = {■^Jp'j + m|, Pj) in the centre-of-mass frame of the hard process, we want to 
preserve the total energy in this frame. 



^' = J2\/'^'+Ph (3-24) 

while we want to keep the sum of spatial momenta vanishing. As the jet parents have 
momenta qj = {Ej, qj) after the showering, we need some way to restore momentum 
conservation in a way that most smoothly preserves the internal properties of each 
jet. 

The simplest way to do this so-called "momentum reshuffling" is to rescale the 
momentum of each jet with a common factor k that is determined implicitly from 



= E V M • (3.25) 

Then, for every jet we determine a Lorentz transformation such that 

q, = {E,,q,)^-^q'^ = iE'^,kp,). (3.26) 

Typically the rescaling parameters k are very close to unity and hence the resulting 
angles and boost parameters in eq. ( |3.26| ) are small. 



4. The cluster hadronization in detail 

The cluster hadronization has two main steps. The first is the cluster formation, 
in which all of the colour connected partons created in the shower are combined to 
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form clusters which are colour singlets. The other step is the decay of these colour 
singlet clusters into hadrons. The new model presented here only changes the second 
stage, how the cluster decays. The process of cluster formation remains identical to 
HERWIG, with the same set of parameters. 



4.1 Cluster formation 



The gluons in the partonic final state are split non-perturbatively into qq pairs. The 
choice of flavour is between the m, d and s flavours. The splitting is done with a 
simple isotropic decay where the gluon is given an effective gluon mass, > 2mq. 



The default value for nig is 0.75 GeV. 



Once we have a state of all on-shell quarks, the colour partners are combined into 



clusters. Owing to the colour-preconfinement property of the parton shower [|T8[, the 
cluster mass distribution is independent of the nature and energy of the hard process 
to a good approximation. This can be seen in Figure ^ for light [uds) quark clusters 
and clusters containing a h quark separately. 



Primary Light Clusters 



Primary fo-Clusters 



0.9 
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0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





Q = 35 GeV 
Q = 91.2 GeV 
Q = 189 GeV 
Q = 500 GeV 
Q= 1000 GeV 
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Figure 2: Primary cluster mass distribution in e^e annihilation at various centre-of-mass 
energies Q for clusters containing only light quarks (left) and a b quark (right). 



The hadronization model in HERWIG and Herwig++ also has a stage where some 
of these clusters are decayed into two new clusters, rather than directly to hadrons. 
This step is called cluster fission. The mass of a cluster is given by = p^, where 
p is the momentum of the cluster. The cluster C is decayed into two new clusters 
Ci, 6*2 if this mass does not satisfy the condition 

^Clp„„ ^ Cl^ax^^""" + (4.1) 

where Clpow and Clmax are parameters of the model and Sc is the sum of the masses 
of the constituent partons which form the cluster. If a cluster does decay into two 
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new clusters, a flavour is drawn from the vacuum. Again this is drawn from the u, d 
and s flavours. The mass of cluster i is drawn from the following distribution 

M, = [(M^ - K + "^3)'') + K + m^Y] , (4.2) 

where rrii is the mass of the constituent quark from the original cluster that is going 
into the new cluster and is the constituent mass of the flavour that was drawn 
from the vacuum. Here P is a parameter of the model and is a random number. 
The value of P is given by PSPLT(l) if parton i is of u,d,s or c flavour and by 
PSPLT(2) for a cluster where parton z is of 6 flavour. The two masses are also 
correlated by the constraint that Mi + M2 < M. If this constraint is violated, a new 
flavour is drawn from the vacuum and two new cluster masses are drawn from the 
distribution ( |4.2| ). The decay kinematics is determined in the rest frame of C, as the 
original constituent quarks continue their movement in the same direction also when 
they are boosted into the rest frame of the new clusters Ci,C2- As all masses are 
given, the momenta of clusters and constituents are determined. 

4.2 Cluster decays 

The last stage of the hadronization is the cluster decays. The problem with the 
original HERWIG cluster decay model can be shown as follows. The probability 



of accepting a decay of a cluster with flavours i,j into hadrons of type a, b is 

P(a,,„6,,,K,j) = PgP{a\i,q)P{b\q,j)Pps{a,b) . (4.3) 

Here Pq is the probability of drawing flavour q from the vacuum and Pps is the 
probability due to phase space. The probabilities of interest are the other two. 
These have the form 

P^a\t,j) = (4.4) 

where Wa is a hadron specific weight, Mij is the maximum weight of all the hadrons 
of flavour i,j and Nij is the number of hadrons of flavour i,j in the model. We 
can see that the probabilities have a dependence on N^j. As described in |T2| this 



causes adverse side effects when new hadrons are added to the model. If we consider 
adding a new hadron of flavour ud, for example, it will generally be heavier than 
those already present. This will suppress the probability of choosing a lighter ud 
meson as A^„j is increased. Therefore properties such as the charged to neutral pion 
ratio are controlled heavily by how many hadrons of a particular flavour are in the 
model. 

To solve this problem a new construction of the probability was created in [|12 



Instead of independently choosing the flavour from the vacuum and then choosing 
the hadrons, this is all combined into one distribution. The weight of one choice is 

W{ai^q, bgj\i,j) = PgWaWbPpsia, b) . (4.5) 
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This gives the probabiUty 

Because Pps is zero for heavy decay modes only accessible modes influence the 
probabilities. Unfortunately, this solution has a new problem in that the ratio of 
mesons to baryons is dictated by the number of available particles. Since there are 
many more mesons then baryons the denominator in eq. ( [4. 61 ) is quite large and the 
total probability of choosing a baryonic decay mode is very small. So though this new 
approach is able to make quantities such as pion ratios independent of the number 
of hadrons in the model it fails to produce the correct amount of baryons. 

The solution implemented in Herwig++ is to treat the baryon sector indepen- 
dently from the meson sector. This is done by re-interpreting the parameter for the 
diquark weight, Pwtdi, to be the parameter that controls the frequency of drawing 
independently from the baryon sector. This is expressed as 

Pb = (4.7) 
Pwtdi + 1 

So there is a probability Pb of choosing only from the baryon sector and a probability 
Pm = 1 — Pb of choosing from the meson sector. The actual choice of hadrons is 
then made according to the probability 

P{ai,q,Ogj\i,j) - — — — — -. (4.8) 

22M/BWiCi,g',dg',j\t,3) 

where the sum over M/B indicates only summing over the flavours that produce 
either mesons or baryons. 

4.3 Hadron decays 

Most of the hadrons created in the cluster hadronization are not stable and need to 
be decayed. At present, the decays in Herwig++ are done in the same way as in 
HERWIG. Most decays are treated as simple n-body isotropic decays. Weak decays 
are done by either free particle V — A matrix elements or bound quark V — A matrix 
elements. Examples of these are e^Vf.VT and K~ e~z7e7r°, respectively. 

Heavy hadrons, such as B mesons, are sometime decayed into partonic states. 
These states are then fed back into the shower and are re-showered and re-hadronized. 
There are two different types of heavy partonic decays. One is a weak decay, for which 
we use the same free or bound V — A matrix elements as for the light mesons. This 
would occur, for example, in the decay 'seed. There are also quarkonium 

decays that have gluons as decay products, for example the decay rjc — > gg. These 
decays are done using the appropriate positronium matrix elements. 
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5. Observables 



5.1 Observables considered 

We considered the following observables in our study. 

Event shape variables: Event shape distributions have been measured to very 
high accuracy at LEP and aim at resolving the properties of the parton shower quite 
deeply. In particular we have chosen the thrust (T), thrust major (M), thrust minor 
(m) and oblateness (O) as they are the most commonly used. In addition we look 
at the C-parameter and D-parameter as they are more sensitive to multijet events. 
We also look at the sphericity and planarity {S,P), which are calculated from the 
quadratic momentum tensor and therefore put more emphasis on 2-jet like events. 
Furthermore we look at the wide and narrow jet broadening measures (-Bmax, -Bmin), 
which are more sensitive to the transverse jet structure. 

Jet multiplicity: The multiplicity of (mini-)jets in e^e~-collisions for different 



values of the jet resolution ?/cut- We use the Durham- or fc^-clustering scheme |19 
throughout the paper for jet observables. To be specific, for a given final state the 
jet measure 

2mm{Ef,E^) 

= — -Q2 — " ^°^%) (^-1) 

is calculated for every particle pair (i, j). The particles with minimal jet measure 
are clustered such that the momentum of the clustered pseudo-particle is the sum 
of the four-momenta of the constituents. The jet multiplicity is then the number of 
pseudo-particles remaining when all Uij > ycut- This inclusive observable has been 
predicted and measured at LEP energies and will test the dynamics of the parton 
shower as well as the interface between parton shower and hadronization. We use 



the Kt Jet-package that implements the above jet-finding algorithm in C++ and 
have written a simple wrapper around it in order to use it with our own particle 
record. 

Jet fractions and 1^: A closer look 'into' the jets is provided by considering the 
rates of jets at a given value of ?/cut in the Durham scheme, Rn = A^„_jet/A^evts for 
n = 2 up to 72 = 6 jets. We also look at the distributions of F„, the ?/cut-values at 
which an n + 1 jet event is merged into an n-jet event in the Durham clustering 
scheme. Here we look at n = 2 up to n = 5. These distributions will not only probe 
the dynamics of the parton shower but also the hadronization model: at the lowest 
values of ?/cut ~ {Qc/QY the dynamics is dominated by the latter. 

Four— jet angles: A very interesting set of observables are the distributions of the 
angles between jets in four-jet events, xbz, ^ksw, ^nr 034, defined for example 



in ref. |21]. These angles are expected to be sensitive to the accuracy of the simulation 



of higher-order matrix elements. 
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Single particle distributions: Another interesting set of observables are the mo- 
mentum distributions of final-state particles with respect to the event orientation. 

is the rapidity distribution with respect to the thrust axis, in and are 
the respective transverse momenta within and out of the event plane, defined by the 
thrust and thrust major axes. Without any reference to the event orientation, we 
look at the distribution of the momentum fraction Xp — 2\p\/Q and at = — log 
which displays better the effects of soft gluon coherence at small Xp. 

Identified particle spectra: We consider the exclusive momentum distributions 
of TT^, K^, p,p and A, A. These are generally expected to be sensitive to the 
hadronization model. In all cases except A, A we can compare with data on the 
momentum distributions from uds, c and b events separately. 

Hadron multiplicities: The charged particle multiplicity distribution and the 
average multiplicities of a wide range of hadron species were taken to test the overall 
fiow of quantum numbers through the different stages of simulation. The improved 
cluster hadronization model can be tested thoroughly against these observables. 

B fragmentation function: The energy fraction of S-hadrons is taken as a test 
for the new parton shower which is specifically designed to improve the description 
of heavy quark observables with respect to the description in HERWIG. 

The above list of observables has proven to be very useful to test different do- 
mains of the available phase space of parameters and has led us to important con- 
clusions for the ongoing development of the code for hadronic collisions. 

5.2 Analysis 

We have booked histograms for all the above distributions in the same bins as the 
experimental data. For a given bin i we then compare the data Di value with 
the Herwig++ Monte Carlo result Mj. Given the data errors 6Di (statistical plus 
systematic) and Monte Carlo errors 5Mi (statistical only), we can calculate a for 
each observable. We keep the statistical error of the Monte Carlo generally smaller 
than the experimental error. In distributions where the normalization is not fixed, 
such as momentum spectra, we allow the normalization of the Monte Carlo to be 
free to minimize x^- The normalization is then tested separately against the average 
multiplicity. In all other cases we normalize histograms to unity. 

As we do not want to put too much emphasis on a single observable or a particular 
region in phase space where the data are very precise, in computing we set the 
relative experimental error in each bin to ma.x{SDi/ Di, 5%). This takes into account 
the fact that the Monte Carlo is only an approximation to QCD and agreement with 
the data within 5% would be entirely satisfactory. The general trend for the preferred 
range of a single parameter was however never altered by this procedure. 
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After normalization the ratio 



is computed for each bin in order to see precisely where the model fails. This ratio 
as well as the relative experimental error (yellow) and the relative contribution of 
each bin to the of an observable is plotted below each histogram. 

5.3 Strategy 

We have taken values for hadron multiplicities into account in the same way as we 
weighted the event shapes. In general the multiplicities of individual particle species 
are sensitive to a completely different set of parameters. The general strategy was 
to start from an initial set of hadronization parameters taken from HERWIG, and 
to aim for a good value for the total number of charged particles with reasonable 
values for the parton shower cutoff parameter b and the maximum cluster mass pa- 
rameter CLMax. Once these were fixed, the hadronization parameters that determine 
the multiplicities of individual particle species were determined. Following this we 
compared this 'preferred' set of parameters with the 'initial' set from HERWIG. The 
resulting parameter set is shown in table |I|. 

6. Results 

We have chosen a wide range of observables in order to test different aspects of the 
model. Event shape variables and multiplicities are considered in order to test the 
dynamical aspects of parton shower and hadronization models, which are closely 
linked at their interface via the parton shower cutoff parameter 5. Ideally, the two 
models should merge smoothly at scales where ~ 1 GeV. All figures shown at the 
end of the paper contain three sets of plots: 

• the actual distribution. The Herwig++ result is plotted as a histogram together 
with the experimental data points (top); 

• the ratio Ri ( ^.2|) together with an error band, showing the relative experimen- 
tal statistical and systematic errors (middle); 

• the relative contribution of each data point to the total of each plot (bottom). 
6.1 Hadron multiplicities 

Table ^ shows the results of the new algorithm in comparison with the old algorithm. 
The column labelled 'Old Model' is the result of using the old algorithm with the new 
shower variables. The column Herwig++ is using the new algorithm with the new 
shower and the last column, labeled Fortran, is using the Fortran HERWIG program 
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(version 6.5). Data are combined and updated from a variety of sources, see ref. [Q. 
We see that, even before systematic tuning, the overall results are better than those 
of HERWIG, with fewer prediction that differ from the data by more than three 
standard deviations (starred). 

We also considered the distribution of the charged particle multiplicity in com- 
parison to OPAL data and find fairly good agreement (fig. |^), although with 



some excess at low multiplicity. 
6.2 Jet multiplicity 

In fig. ^ we show the average number of jets {rtjets) at the Z°-pole, as a function of 
the Durham jet resolution Ucnti for various values of the cutoff parameter 5. At the 
parton level (top left) the jet multiplicity varies a lot as we go to small values of ycnti 
saturating at the number of partons that are present in a single event. The order 
of magnitude of the visible saturation scales is characterized for each flavour by the 
different cutoff values Qg as y^ax = Q\IQ'^ (cf- ( p.l6|) ). During hadronization, low 
parton multiplicities lead to large mass clusters which tend to decay into low mass 
clusters below a given cutoff mass, which is fixed to its default value throughout 
the current section. Fig. ^ (top right) shows that the hadronization compensates for 
lower partonic multiplicities, giving a result insensitive to 5 at the hadron level. In 
other words, we have a smooth interface between perturbative and non-perturbative 
dynamics between the lower end of the parton shower on one side and the cluster 
hadronization model on the other side. On the hadron level we describe LEP data 



from OPAL |2J well. 

In order to test the sensitivity of our model against the variation of the cm. 
energy, we calculate the jet multiplicities at PETRA and LEP II energies as well 



(fig. ^, bottom). The comparison to JADE |]23[ and OPAL data shows a good 
agreement. In all runs we applied the same cutoffs on the energy of the partonic 
subsystem as did the experiments. 

The additional curves in fig. ^ show predictions for the jet multiplicity from 



the resummation of leading logarithms. Note that the parameter Aqcd in the re- 
summed calculation is not Aj^. We see that for the value Aqcd = 500 MeV there is 
good agreement with the data and the Herwig++ result throughout the perturbative 
region, |/cut > 10""^. 

6.3 Jet rates and Yn distributions 

Another set of observables that is known to be well-described at LEP energies are the 
fractions of n-jet events at a given ?/cut in the Durham scheme. In fig. ^ we compare 
the results from Herwig++ with LEP data from |Q and find good agreement. On 



the hadron level these predictions are not very sensitive to the cutoff parameter 5. 
The 5-jet distribution is not shown and R% is the rate of > 6-jet events. 
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The Durham Yn distributions in fig. ^ are histograms of those |/cut ^values at which 
an n + 1-jet is merged into an n-jet event in the Durham jet clustering scheme. We 
may say that they resolve more internal structure of the jets than the ra-jet rates. 
Still, the agreement between model and data is quite good, although there is a 
tendency (which was also present in HERWIG) to exceed the data at low Yn. 



6.4 Event shapes 

In order to test the dynamics of the parton shower in Herwig++ in more detail we 
consider a set of commonly used event shape variables. Not only the coUinear region 
of the parton shower is probed in greater detail but also the regions of phase space 
which are vetoed as matrix element corrections. We compare all results to DELPHI 
data |2^. 

In fig. I^we show the distribution of thrust, thrust-major and thrust-minor. These 
variables are all obtained from a linear momentum tensor. The thrust distribution 
is shown with and without matrix element corrections switched on. The prediction 
without matrix element corrections is very much better than that of HERWIG, 
owing to the improved shower algorithm. It is interesting that the matrix element 
corrections seem to generate almost too much transverse structure, leading to event 
shapes that are less two-jet-like. On the other hand, there is also a slight excess of 
events close to the two-jet limit. 

It is remarkable how well distributions like C and D parameter (fig. ^ which are 
sensitive to three- and four-jet-like events are described by our model even though 
we are limited to three jet matrix elements plus showers. Here again we have in fact 
a small excess at high values. 

We show also in fig. ^ the distributions of sphericity and planarity, which are 
obtained from a quadratic momentum tensor and therefore put more emphasis on 
high momenta. As was the case for the thrust-related distributions, we tend to have 
slightly wide events. In addition we consider the jet broadening measures -Bmax and 
-Bmin and the hemisphere jet masses (fig. ||). In all cases the agreement between 
model and data is good. 



6.5 Four jet angles 

We show the four-jet angles in fig. They are considered only for events where we 
have a four-jet event at i/cut = 0.008. Despite the fact that we do not have any match- 
ing to higher order matrix elements, as was proposed in p8| and implemented in [0], 



the agreement between model and data [29| is remarkably good. Even though we ex- 
pected the implementation of hard and soft matrix element corrections in Herwig++ 
to be most important for the description of these observables, we did not find very 
significant differences with or without the application of matrix element corrections. 
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6.6 Single particle distributions 



In fig. |IT] we show single charged particle distributions within the event, oriented 
along the thrust axis. The transverse momentum within the event plane is 
shown with and without matrix element corrections. In contrast to the thrust distri- 
bution we find that the matrix element corrections actually improve the distribution. 
Furthermore, and the rapidity along the thrust axis are rather well described. 

We do not show the analogous momentum distributions with respect to the sphericity 
axis which have similar features. 

We consider the distribution of scaled momentum Xp = 2\p\/Q of charged par- 
ticles in fig. |T2|. In addition to the full distribution we also consider the results from 
light (uds), c and b events.'^ In all cases we compare with data from SLD The 



charged particle distribution is well described in all four cases, in fact somewhat 
better for heavy primary quarks. 

6.7 Identified hadron spectra 

As in the case of all charged particles we can compare identified particle spectra from 
events of different flavour to SLD data [Q. Data for tt^ (not shown, being almost 
equivalent to all charged particles), K"^ and {p,p) are available. In flg. |13| we see 
the data for {p, p) spectra from events of different flavour. For large values of Xp we 
clearly overshoot the data in light flavoured events. This is somewhat compensated 
by the heavy quark events which in turn seem to prefer lower values of Xp. We believe 
that this feature is related to the hadronization, being similar to but smaller than 
that seen in HERWIG. 

Fig. |14| shows distributions for and A, A. Both are rather better described 
than the proton spectra but the distribution of A, A tends to have a similar, though 



smaller, 'bump' in comparison to data from ALEPH 
6.8 B fragmentation function 



In flg. 113] we consider the B hadron fragmentation function in comparison to data 



from SLD WM. We have also considered data from ALEPH P3| (not shown). We 



can describe the data quite well without any additional tuning of the hadronization 



model to this data. The parton shower formulation in terms of the new variables |]10 
and taking quark masses in the splitting functions into account clearly improves the 
description of heavy quark events. 

6.9 Overall results 

In table ^ we show a list of values for all observables that were studied during 
our analysis, including those not shown in the plots. The most sensitive parameters 
were the cutoff value 6 and the use of (hard plus soft) matrix element corrections. 

^The flavour of the quark-antiquark produced in tlie initial hard process. 
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The table shows three values of 5: our preferred value of 5 = 2.3 GeV as well as the 
lowest and highest values that we considered. 

The results should be interpreted with care. The overall trend suggests that we 
should prefer a large cutoff scale. However, we have just averaged over all possible 
observables. Taking a closer look, we may want to weight different observables in a 
different way. 

In more detail, the general trend is the following: event shapes, jet rates and 
differential jet rates prefer a low cutoff. The single particle distributions along the 
thrust and sphericity axes prefer a small cutoff value. The ynm distributions prefer 
either a high or a low cutoff value. The spectra of identified particles tend to prefer the 
high cutoff value with some exceptions for light quark events. The B fragmentation 
function clearly prefers the intermediate value. 

In addition, as indicated in sec. ^]3|, we found that the measured yields of iden- 
tified particles clearly prefer the value 6 = 2.3 GeV. 



7. Conclusions 

We have achieved a complete event generator for e^e~ annihilation into hadrons. 
The main physics features, in comparison to the previous versions of HERWIG, are 
an improved parton shower, capable of properly describing the perturbative splitting 
of heavy quarks, and an improved cluster hadronization model. 

We have tested our model against a wide range of data from e~^e~ colliders and 
are able to give a good general description of the data. 

For many observables the description of the data has been improved with respect 
to HERWIG. The new parton shower has a number of remarkable features. The need 
for matrix element corrections has decreased. The main reason for this is the use 
of improved splitting functions, which give a far better approximation of the matrix 
elements in the region of collinear gluon emissions. We can describe observables 
involving light or heavy quark splitting with a unique set of parameters. The new 
hadronization model also improves the description of identified particle spectra and 
multiplicities. 

The detailed analysis of our results leaves us with a recommendation: the set 
of parameters that is shown in table |I|. This set of parameters is understood as a 
weighted compromise in order give a good overall description of the data we have 
considered so far. We did not aim at a complete tuning of the model, but rather 
wanted to study its ability to describe the broad features of the data, which turned 
out to be very successful. 

Future work on the program will extend the parton shower to initial state radi- 
ation and include a model for the soft underlying event in hadron-hadron collisions, 
aiming at a complete event generator for the simulation of Tevatron and LHC events. 
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Parameter 


Preferred 


Initial 


as{Mz) 


0.118 


0.114 


S/GeV 


2.3 


— 


nig/GeV 


0.750 


— 


Qmin/GeV in ^^(Qinin) 


0.631 


— 


CW/GeV 


3.2 


3.35 


^-■-pow 


2.0 


— 


PSpltl 


1 


— 


PSplt2 


0.33 


— 


BlLim 


0.0 




ClDirl 


1 




ClDir2 


1 




ClSmrl 


0.40 




ClSmr2 


0.0 




Pwtd 


1.0 




Pwt^j 


1.0 




Pwt^ 


0.85 


1.0 


Pwtc 


1.0 




Pwtfe 


1.0 






0.55 


1.0 


Singlet Weight 


1.0 




Decuplet Weight 


0.7 


1.0 



Table 1: The preferred parameters for Herwig++. The first group are shower parameters, 
the second are all of the hadronization parameters. In the third column we show initial 
values of our study, taken from HERWIG, where these differ from the preferred values. 
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Particle 


Experiment 


Measured 


Old Model 


Herwig++ 


Fortran 


All Charged 


M,A,D,L,0 


20.924 ± 0.117 


20.22* 


20.814 


20.532* 


7 


A,0 


21.27 ± 0.6 


23.03 


22.67 


20.74 




A,D,L,0 


9.59 ± 0.33 


10.27 


10.08 


9.88 


/9(770)0 


A,D 


1.295 ± 0.125 


1.235 


1.316 


1.07 


7r± 


A,0 


17.04 ± 0.25 


16.30 


16.95 


16.74 


/3(770)± 





2.4 ± 0.43 


1.99 


2.14 


2.06 


?? 


A,L,0 


0.956 ± 0.049 


0.886 


0.893 


0.669* 


a;(782) 


A,L,0 


1.083 ± 0.088 


0.859 


0.916 


1.044 


r/'(958) 


A,L,0 


0.152 ± 0.03 


0.13 


0.136 


0.106 




S,A,D,L,0 


2.027 ± 0.025 


2.121* 


2.062 


2.026 


K*(892)0 


A,D,0 


0.761 ± 0.032 


0.667 


0.681 


0.583* 


E:*(1430)° 


D,0 


0.106 ± 0.06 


0.065 


0.079 


0.072 




A,D,0 


2.319 ± 0.079 


2.335 


2.286 


2.250 




A,D,0 


0.731 ± 0.058 


0.637 


0.657 


0.578 


(/>(1020) 


A,D,0 


0.097 ± 0.007 


0.107 


0.114 


0.134* 




A,D,0 


0.991 ± 0.054 


0.981 


0.947 


1.027 


A++ 


D,0 


0.088 ± 0.034 


0.185 


0.092 


0.209* 


s- 





0.083 ± 0.011 


0.063 


0.071 


0.071 


A 


A,D,L,0 


0.373 ± 0.008 


0.325* 


0.384 


0.347* 


EO 


A,D,0 


0.074 ± 0.009 


0.078 


0.091 


0.063 


E+ 





0.099 ± 0.015 


0.067 


0.077 


0.088 


S(1385)=^ 


A,D,0 


0.0471 ± 0.0046 


0.057 


0.0312* 


0.061* 




A,D,0 


0.0262 ± 0.001 


0.024 


0.0286 


0.029 


H(1530)0 


A,D,0 


0.0058 ± 0.001 


0.026* 


0.0288* 


0.009* 


0- 


A,D,0 


0.00125 ± 0.00024 


0.001 


0.00144 


0.0009 


/2(1270) 


D,L,0 


0.168 ± 0.021 


0.113 


0.150 


0.173 


/^(1525) 


D 


0.02 ± 0.008 


0.003 


0.012 


0.012 




A,D,0 


0.184 ± 0.018 


0.322* 


0.319* 


0.283* 


L»*(2010)=^ 


A,D,0 


0.182 ± 0.009 


0.168 


0.180 


0.151* 


DO 


A,D,0 


0.473 ± 0.026 


0.625* 


0.570* 


0.501 




A,0 


0.129 ± 0.013 


0.218* 


0.195* 


0.127 







0.096 ± 0.046 


0.082 


0.066 


0.043 




A,D,L,0 


0.00544 ± 0.00029 


0.006 


0.00361* 


0.002* 




D,0 


0.077 ± 0.016 


0.006* 


0.023* 


0.001* 


^''(3685) 


D,L,0 


0.00229 ± 0.00041 


0.001* 


0.00178 


0.0008* 



Table 2: Multiplicities per event at 91.2 GeV. We show results from Herwig++ with the 
implementation of the old cluster hadronization model (Old Model) and the new model 
(Herwig++), and from HERWIG 6.5 shower and hadronization (Fortran). Parameter 
values used are given in table ||. Experiments are Aleph(A), Delphi(D), L3(L), Opal(O), 
Mk2(M) and SLD(S). The * indicates a prediction that differs from the measured value by 
more than three standard deviations. 



- 25 - 



1 I r 

Herwig++ 1.0 

6 = 1.7 GcV 

5 = 2.3GeV 

6 = 3.2 GcV 

OPAL 99 




10 20 30 40 50 



Figure 3: The distribution of the charged particle multipUcity. The three panels here and 
in figs. 5-15 are explained in sec. 6. 



- 26 - 



Partons 



Hadrons 





lQ-5 IQ~4 iqO iQ-b ]^Q-:i ;^Q-2 ;^g( 



R2 vs ;y„,t R-j vs j/cut 




10-= 10-** 10-'* 10-^ 10-1 10" 10-= lo-'' iQ-^ 10-^ 10-1 10" 



Figure 5: Jet rates in tlie Duriiam algoritlim for different values of the cutoff 6. 
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Figure 6: Durham Yn distributions for different values of the cutoff 6. 
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Figure 7: Thrust without (top left) and with (top right) matrix element corrections 
switched on, thrust major and thrust minor (bottom). 
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Figure 8: Sphericity, planarity, C parameter and D parameter distributions. 
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Figure 9: The wide and narrow jet broadening measures Smax and i?min and the high 
and low hemisphere masses. 
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Figure 10: Four jet angle distributions. The points are from preliminary DELPHI data. 
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Figure 11: Momentum distributions of charged particles with respect to the thrust axis, 
Pi. in (with and without matrix element corrections), and y'^ . 



-34- 




Figure 12: The scaled momentum distribution Xp of charged particles for all events as 
well as for uds, c and h events separately. 
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Figure 13: The scaled momentum distribution Xp of protons, shown separately for all 
events as well as for uds, c and b events. 
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Figure 14: Distribution of scaled kaon momentum and A, A momentum. 
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Figure 15: The i?-hadron fragmentation function for different values of the cutoff 6. 
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ME corrections off ME corrections on 



Observable Ref. 


.5 = 1.7 GcV 


<5 = 2.3 GeV 


(5 = 3.2 GeV 


<5 = 1.7 GeV 


<5 = 2.3 GcV 


<5 = 3.2 GeV 


1 - T 


27 




44.65 


33.15 


22.29 


72.80 


45.57 


26.34 


M 


27 




246.25 


273.42 


198.37 


275.80 


274.43 


186.34 


m 


27 




150.74 


157.91 


137.43 


174.29 


163.02 


129.10 


O 


27 




7.41 


5.58 


5.14 


22.24 


19.33 


13.34 


s 


27 




4.42 


3.50 


4.07 


24.70 


14.10 


8.50 


p 


27 




4.48 


5.63 


6.54 


10.69 


7.40 


5.62 


A 


27 




19.52 


10.80 


7.17 


44.83 


20.26 


11.33 


C 


27 




66.86 


59.26 


39.56 


81.41 


67.44 


43.08 


D 


27 




84.23 


29.30 


12.36 


161.90 


60.42 


26.92 


A-^high 


27 




25.78 


18.88 


12.38 


38.43 


25.69 


11.52 


Mlow 


27 




15.25 


5.37 


2.50 


31.53 


10.42 


5.00 


Mdift 


27 




7.28 


5.27 


7.25 


18.32 


12.17 


4.61 


^max 


27 




54.48 


50.29 


38.91 


59.61 


49.92 


33.23 


-^min 


27 




53.25 


55.72 


53.18 


64.52 


58.08 


50.64 


^sum 


27 




102.29 


97.35 


74.60 


121.86 


103.10 


70.98 


Bdiff 


27 




8.28 


5.42 


4.70 


18.39 


13.64 


6.09 




27 




2.48 


3.11 


11.52 


3.39 


1.70 


4.26 


T ' 
Pj_ out 


27 
27 




0.25 
34.52 


3.28 
60.55 


21.65 
66.05 


0.80 
34.94 


1.70 
53.81 


16.06 
59.07 


flin 


27 




2.53 


3.19 


11.76 


2.32 


1.39 


4.30 




27 
27 




0.37 
9.04 


3.77 
17.49 


22.64 
24.85 


0.90 
7.78 


2.01 
14.72 


16.78 
21.94 




27 




9.37 


3.54 


3.76 


25.56 


11.27 


5.25 


3 


27 




25.85 


6.33 


2.14 


47.11 


15.31 


5.42 


i^F 

4 


27 




43.90 


10.47 


2.69 


78.82 


23.26 


7.11 




24 




8.75 


6.11 


5.36 


12.35 


8.65 


6.40 


7/34 


24 




10.20 


9.65 


9.07 


11.46 


10.02 


8.81 




24 




15.53 


14.40 


11.78 


17.74 


15.57 


11.75 


ysG 


24 




16.02 


17.77 


15.13 


15.50 


17.51 


14.32 


(A^jcts) 


24 




12.84 


3.30 


0.62 


28.29 


12.80 


5.95 




24 




9.75 


6.56 


6.18 


19.84 


13.45 


9.59 




24 




10.46 


8.51 


9.36 


23.49 


15.86 


11.95 


i?4 


24 




13.47 


10.95 


10.36 


15.26 


12.42 


10.22 




M 




25.53 


24.98 


23.43 


28.09 


26.35 


22.30 


-Re 


24 




10.37 


1.74 


0.67 


18.38 


4.33 


1.47 


^^'^xAdZj ) 


2g 




2.90 


1.10 


0.48 


2.48 


1.05 


0.53 


COsf '^'t^qw^ 


29 




2.30 


2.06 


2.56 


1.22 


1.50 


1.64 


V JN H' 
COs(q34) 


2S 
2S 




7.68 
1.41 


5.06 
1.57 


2.72 
1.71 


8.66 
0.60 


6.22 
0.64 


3.57 
0.76 




23 




21.86 


25.71 


12.90 


19.81 


22.84 


12.97 


3;p{ch)[all] 


3C 




5.32 


5.65 


3.49 


4.75 


4.10 


3.02 


a:p(ch)[uds] 


3C 




15.72 


8.50 


6.13 


12.63 


6.69 


5.86 


(ch) [c] 


3C 




3.95 


2.29 


2.17 


2.96 


1.76 


2.73 


a^nfch) [61 


3C 




35.05 


3.23 


1.79 


35.79 


2.49 


1.22 


z„('7r±Malll 

V V / L J 


3C 




8.29 


9.27 


6.18 


7.21 


7.50 


5.51 




3C 




28.30 


15.92 


10.47 


24.05 


13.29 


9.46 


f \ J L J 


3C 




4.65 


2.99 


1.38 


3.67 


2.28 


1.62 


xJ-K^)\b] 

/L J 


3C 




49.13 


3.14 


1.56 


49.44 


3.57 


2.02 


Xp{/s:±)[all] 


3C 




4.99 


2.02 


15.38 


3.67 


2.88 


17.37 


xp(_fs:±)[uds] 


3C 




6.46 


17.05 


36.45 


6.83 


19.36 


38.79 


x^(K±)[b\ 


3C 




21.01 


2.20 


3.35 


18.16 


1.75 


4.14 


3C 




8.56 


7.14 


4.34 


7.63 


5.84 


4.97 


2:p(p>p)[all] 


3C 




143.34 


98.19 


42.90 


140.48 


87.08 


36.23 




3C 




145.35 


102.51 


52.78 


139.85 


91.07 


45.10 


^p(P.P)[c] 


3C 




2.26 


2.41 


2.86 


2.34 


2.48 


2.85 




3C 




11.26 


13.71 


8.12 


11.47 


13.54 


8.31 


p{A,A) 


31 




58.02 


28.52 


9.47 


55.27 


25.50 


7.86 




32 




8.93 


0.92 


8.16 


9.44 


1.39 


9.92 




33 




15.40 


1.75 


7.35 


15.76 


2.01 


8.21 


(x">/bin 






32.75 


25.84 


20.93 


40.69 


28.41 


20.56 



Table 3: x^/bin for all observables we studied and a relevant subset of parameters. 
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